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C^l ■ Abstract. We consider partial sum rules for the homogeneous limit of the solution of 

(— 1 ■ the q-deformed Knizhnik-Zamolodchikov equation with reflecting boundaries in the Dyck 

i path representation. We show that these partial sums arise in a solution of the discrete 

Hirota equation, and prove that they are the generating functions of r 2 -weighted punc- 
£ — ' tured cyclically symmetric transpose complement plane partitions where r = — (q + q^ 1 ). 

In the cases of no or minimal punctures, we prove that these generating functions coincide 
with T 2 -enumerations of vertically symmetric alternating sign matrices and modifications 
' thereof. 
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1. Introduction 

The discrete Hirota equation [29] is a difference equation for a function / in three 
discrete variables: 

f(n,i,j)f{n - 2,i,j) = /(n - l,i - 1, j)f(n -l,i+ 

+ r 2 f(n - - l)f(n - + 1). (1.1) 

The factor r 2 can be absorbed into / by rescaling f(n,i,j) — > r _J ^ +1 )/(n, i, j). Equa- 
tion (|l.ip is also known as the octahedron recurrence as the variables n,i,j are natural 
coordinates for the vertices of an octahedron, see e.g. [26]. 

In order to fix a particular solution to (II. ip we need to specify boundary conditions. 
An interesting choice is the following. Let A = (ajj)i<jj< n be an n x n matrix and set 

/((Ui)=l, f(l,i,j)= aij l<i,j<n. (1.2) 

With these boundary conditions, the discrete Hirota equation defines the r 2 -deformation 
of the determinant studied by Robbins and Rumsey [25]. In particular, the r 2 -determinant 
of A is defined by 

L4| T 2 = /(n,l,l). (1.3) 

For r = i the r 2 -determinant reduces to the ordinary determinant which can be expanded 
as a sum over permutation matrices. In their famous work [25], Robbins and Rumsey 
showed that solutions to (jl.ip with the boundary condition (jl.2p can be written as a sum 
over alternating sign matrices (ASMs), see also [2]. 

It is well known that ASMs are equinumerous to totally symmetric self-complementary 
plane partitions (TSSCPPs). Surprisingly, generating functions of T 2 -enumerations of 
TSSCPPs and other symmetry classes of plane partitions, as studied by Robbins |24j . 
appear as normalisations for homogeneous solutions of the (/-deformed Knizhnik-Zamolod- 
chikov (qKZ) equation recently obtained by Di Francesco and Zinn-Justin (3(9]. This result 
will be generalised below, when we consider a class of punctured cyclically symmetric 
transpose complement plane partitions (PCSTCPPs), whose weighted enumerations also 
arise in the qKZ equation. 

It was already observed by Robbins and Kuperberg \2i\ [13] that the r 2 -enumerations of 
CSTCPPs (without puncture) are closely related to r 2 enumerations of vertically symmet- 
ric alternating-sign matrices (VSASMs) and other symmetry classes of ASMs. A precise 
statement will be proved below. Closing the circle, these enumerations comprise a partic- 
ular solution of (II. ip . albeit with a different boundary condition than (II. 2|) . We hope that 
this paper will be a further step in resolving the Razumov-Stroganov conjectures and the 
discovery of a bijection between ASMs and TSSCPPs. 

Throughout the following we will use the notation [x] q for the usual q- number 




The notation [x] will always refer to base q. 
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2. (/-DEFORMED KNIZHNIK-ZAMOLODCHIKOV EQUATION 

Definition 1. The Temperley Lieb algebra of type Al, denoted by T^ L (q), is the unital 
algebra defined in terms of generators e,, i = 1, . . . , L — 1, satisfying the relations 

ef = -[2]eu eiej = e^e, Vi, j : \i - j\ > 1, 

e-iei±\ei = e«. (2.1) 

Every representation of the Temperley-Lieb algebra defines a solvable lattice model 
through the definition of so-called R-matrices. An R-matrix is a representation of the 
following Baxterised element of the Temperley-Lieb algebra: 

_ [l-u]- [u] eg . . 

" + ' (2 - 2) 

where u € C \ {— 1} is the spectral parameter. 

Definition 2. A Dyck path a = (ao,«i, • • • ,oll) is a sequence of integers a% (heights) 
such that aj+i = on ± 1, a, > and ao = ckl = 0. 

The Temperley-Lieb algebra has a known action on Dyck paths, which is well-documented 
in the literature, see for example |16} [5] and references therein. The span of Dyck paths 
forms a module of T^ L (q), and we will denote its basis elements by \a), where a runs over 
the set T>l of Dyck paths of length L. Let us now consider a linear combination I*) of 
states | a) with coefficients ip a taking values in the ring of formal series in L variables 
q ±x % i = l,2,...,L: 

\V(x 1 ,...,x L )} = i/j a (xi,...,x L )\a). 
aev L 

The g-deformed Knizhnik-Zamolodchikov equation on a segment (with reflecting bound- 
aries) is a system of finite difference equations on the vector j*). This equation reads (see 

HI) 

Ri(Xi -X i+1 )\iZ>) =7Ti|¥), Vt = l,...,L-l, 

|*)=7T |*>, (2.3) 
|*> = 7TL|tf), 

where R{ are the Baxterised elements of the Temperley-Lieb algebra. The operators 
Ri( 

xi — Xj+i) act on states |q), whereas the operators 7Tj permute or reflect arguments of 
the coefficient functions: 

T^i^cJy- • • i Xi, , • • •) — 1pai\- • • i j %i , • • • ) > 

7T0lp a (xi, . . .) = 1p a (-Xi, . . .), (2.4) 
7T L -0a(- • • , XL) = tpa(- ■ ■ , ~A - X L ). (2.5) 

The shift A S C is a parameter related to the level of the qKZ equation, see [TP] . 

In [7] polynomial solutions of the qKZ equation were studied in the limit where X{ — > 0. 
Interestingly, in this limit the coefficients ip a turn out to be polynomials with positive 
coefficients in r 2 , where r = —[2] = — q — g" 1 . Furthermore, based on observations on 
explicit solutions for small values of L, an intriguing connection between (generalised) sum 
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rules and the enumeration of weighted CSTCPPs was conjectured. This was subsequently 
proved using multi-integral formulae [9]. Here we shall generalise this result. 

3. Explicit solutions 

In [7J explicit solutions of the qKZ equation in the limit X{ — > were obtained for L < 8. 
In order to illustrate our results we shall here list the first few solutions. We will write 
shorthand ijj a for the limit Xj — > of ip a (xi, . . . ,xl). The complete solution is determined 
up to an overall normalisation, and we will choose 

^ n = r LV2j(LV2j-i)/2 

for the coefficient corresponding to the maximal Dyck path 0, E T>l '■ = min{i, L + €l — 
i}, ei = L mod 2. Together with the solution we list certain powers r Ca ' n whose meaning 
will become clear below. 



J — A 


cx 






XX 


l + T 2 1 

r r ±l 


L = 5 


a 


1p a r ±C«,2 r ±C a ,l 




/\/\/\ 

^X> 

xx> 


r 2 (2 + r 2 ) 1 

r(2 + r 2 ) X 1 

1 + 2r 2 t ±2 1 


L = 6 


a 


1p a T ±C a ,2 r ±C a ,l 




/xXX 
XXX 


1 + 5r 2 + 4r 4 + t 6 1 

r(2 + 2r 2 + r 4 ) r ±x 

r(l+3r 2 + r 4 ) X 1 

2r 2 (l + r 2 ) r ±2 1 
r 3 r ±i r ±i 



Table 1. Explicit solutions of the gKZ equation in the homogeneous limit 



An immediate observation about these solutions was already noted in [7]: 

Conjecture 1. The components ip a (xi, . . . , xl) of the polynomial solution of the qKZ 
equation of type A in the limit Xi — > 0, i = 1, . . . , L are, up to an overall factor which is a 
power of t, polynomials in t 2 with positive integer coefficients. Here r = — [2] . 
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Polynomiality and integrality of the coefficients is proved in [9]. 

Based on these explicit solutions, and solutions obtained for L = 9 and L = 10 in [5], 
we discovered underlying discrete bilinear relations. In order to uncover these relations 
we have to introduce certain partial sums over the components ip a of the solution. Let us 
first define the paths £l(L,p) G T>l whose local minima lie at height p, where 

p = L(£-i)/2J -P, p = 0,...,L(£-l)/2j. 

Figure [1] illustrates the path 0(12,3). We further define the subset T>l, p of Dyck paths 
of length L which lie above 0(L,p), i.e. whose local minima lie on or above height p. 
Formally, this subset is described as 

Dl, p = {a G T>l\ cti > Qi(L,p) = mm(Qi,p)} , 

where Oj are integer heights of the maximal Dyck path = 0(L, 0). 



A 

p=3 




L = 12 



Figure 1. The path 0(12,3) G £>i 2)3 . 

To each Dyck path we associate an integer c atP . Let a = (ao, cti, . . . , ol) G T^L,p be a 
Dyck path of length L whose minima lie on or above height p. Then c a p is defined as the 
signed sum of boxes between a and 0(L,p), where the boxes at height p + h are assigned 
(— l)^ 1 for h > 1. An example is given in Figure [21 and an explicit expression for Cq, p is 
given by 

Ca, P = X^" 1 )^ " ^( L >P))- 

i=2 

In the next section we will consider certain properties of the partial weighted sums 

S±(L,p)= r ±c ^i> a . (3.2) 

The partial sums corresponding to the solutions in Table [T] are given in Table 

4. Plane partitions and the discrete Hirota equation 

It was observed in [21] that for r = 1 (q = e 2 "" 1 / 3 ), the partial sums in Table [2] satisfy 
a discrete bilinear relation called Pascal's hexagon, or the discrete Boussinesq equation 
which, in turn, is a two dimensional reduction of the discrete Hirota equation (see [29] and 
references therein). Here we generalise this result to arbitrary r, and show that the partial 
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£2(12,4) 



L=12 



Figure 2. Definition of the number c a)P as the signed sum of boxes be- 
tween the a and the path $7(12,4). In this figure L = 12 and p = 4 and 



Co, 4 - 


= 4- 


-3 + 


1 = 2. 






L = 4 


£_ 


(4,0) 


= r 


£+(4,0) 


= r 




£_ 


(4,1) 


= 2 + r 2 


£+(4,1) 


= 1 + 2r 2 


L = 5 


£_ 


(5,0) 


= r 


5+ (5,0) 


= T 




£_ 


(5,1) 


= 2(l+r 2 ) 


£+(5,1) 


= 1 + 3r 2 




£_ 


(5,2) 


= t- 2 (1 + 5t 2 + 4t 4 + t 6 ) 


£+(5,2) 


= t 2 (6 + 5t 2 ) 


L = 6 


£_ 


(6,0) 


= r 3 


£+(6,0) 


= r 3 




£_ 


(6,1) 


= r 2 (3 + 2r 2 ) 


£+(6,1) 


= r 2 (2 + 3r 2 ) 




£_ 


(6,2) 


= 6 + 13r 2 + 6r 4 + r 6 


£+(6,2) 


= 1 + 8r 2 + 12r 4 + 5r 6 



Table 2. Partial sums corresponding to the solutions in Table [TJ 



sums satisfy the discrete Hirota equation. Based on experimental data for solutions up to 
L = 10 [5], we were led to introduce the following polynomials: 

m + L-2p-k\ _ 2(2m -r) N 




T ^=J$ SP {-L{ r - t ){ \ m -r r J' <41) 

We then have the following result for the functional form of the partial sums for arbitrary 
system size L: 

Proposition 1. 

S + (L,p) = T^T(L,p,[L/2\ -p), (4.2) 
S-(L,p)=T^T(L,p,[L/2] -p+1), (4.3) 



4 The idea to consider such type determinants comes from the observation that some of the polyno- 
mials given in Table [5] coincide with the conjectural generating functions of VSASMs and VHSASMs 
T n (r 2 , n) 1,1=0,1 (see [24], Table 4.3). The latter functions are given by similar determinants. 
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where 



1 



(LX/2j(LL/2j-l)-p(p + l)). 



VL,p 



2 



This proposition will be proved using multiple integral equations in Section [H 

4.1. Punctured symmetric plane partitions. Using the standard presentation of the 
determinant of p x p matrix A in terms of the minors of p x 2p matrices B and C: 
A = BC l , the polynomials T{L,p, k) can alternatively be written as 



By the use of Lindstrom-Gessel-Viennot formula [H] the last expression can be interpreted 
as the generating function of two sets of paths with the same end-point, as in Figure [3l 
One set of paths starts at positions {£,£ + k — 1), {£ = 1, . . . , p) and has diagonal NW-SE 
and vertical steps; the other set starts at {£ — L + 2p + k, — I — L + 2p + k) with diagonal 
NE-SW and vertical steps. Each vertical step below the horizontal line (green online) is 
assigned a weight r 2 . All other steps are assigned weight 1. Up to an overall factor r( fc-1 ) p 
this is also the generating function of pairs of paths where all vertical steps are assigned 
a weight r. 



Figure 3. An example of two sets of paths with common endpoints enu- 
merated by T(L,p, k) for L = 13, p = 4 and k = 3. 

Consider stacking a large cube with smaller cubes, starting in a corner of the large cube, 
in such a way that the smaller cubes can only be stacked on top or next to each other. Such 
stackings are called plane partitions, which are well known to be equivalent to rhombus 
tilings of a hexagon. The paths in Figure [3] arise naturally when one considers cyclically 
symmetric transpose complement plane partitions (CSTCPPs) and generalisations thereof, 
as was done by Ciucu and Krattenthaler [3]. CSTCPPs correspond to rhombus tilings 




r= 1 2 3 4 
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of a hexagon which are invariant under rotations over 2ir/3 as well as under reflections 
across a symmetry axis not passing through the hexagon corners. In [3] the more general 
problem was considered of cyclically symmetric transpose complement tilings of a hexagon 
with a triangular hole, as in Figure HI We will call such tilings punctured CSTCPPs, or 
PCSTCPPs, the size of the triangular puncture being determined by the difference of the 
lengths of the sides of the hexagon. 

If we weight PCSTCPPs by assigning a weight r 2 to each vertical step below the bi- 
secting line (green online) in the South-East region of Figure [H the weighted enumeration 
of tilings of the fundamental domain of PCSTCPPs is equivalent to the weighted enumer- 
ation of paths in Figure We have sketched the bijection in Figure HI full details may 
be found in [3] where the enumeration (r = 1) is considered, which is shown to factorise 
completely, see also f)5. 1 1) . For general r but p restricted to p = \_(L — l)/2j PCSTCPPs 
were considered in the context of the qKZ equation by Di Francesco [7]. 




Figure 4. A punctured cyclically symmetric transpose complement plane 
partition (PCSTCPP) and its fundamental region. The linear size of the 
puncture is given by the difference of the lengths of the sides of the hexagon. 
The position of the bisecting line (green online) determines a particular 
weighting of the PCSTCPP: vertical steps below this line are assigned a 
weight t 2 . 

Interestingly, there is another way to interpret the paths in Figure [3] in terms of punc- 
tured cyclically symmetric plane partitions. When we remove a central cubic region from 
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the large cube, and consider special sixfold rotational symmetric rhombus tilings of the 
corresponding punctured hexagon by defining a fundamental region as in the upper right 
corner of Figure [5j we obtain a subset of punctured cyclically symmetric self complement 
plane partitions (PCSSCPPs). 




L=12 



Figure 5. Fundamental region for a subclass of punctured cyclically sym- 
metric self complement plane partitions (PCSSCPPs). The position of the 
bisecting line (green online) determines a particular weighting of a 
PCSSCPP. 

PCSSCPPs are enumerated by giving a special weight to the non-intersecting lattice 
paths in Figure O see [3]. Here we will not discuss this, but will only consider those 
PCSSCPPs as defined in Figure [6l The paths in this picture correspond to those in 
Figure [3l In the East region of Figure [6l diagonal steps below the bisecting line (green 
online) are assigned a weight r 2 . 

4.2. Discrete Hirota equation. We end this section with the following important ob- 
servation for which we do not know a good interpretation. The polynomials T(L,p,k) 
defined in (|4.ip satisfy a discrete bilinear relation: 

Proposition 2. The polynomials (|4.ip satisfy the recurrence 

T(L,p, k)T(L - 2,p - 2, k + 2) = T(L — l,p — 2,k + 2)T{L - l,p, k) 

+ r 2 T(L - 2,p - 1, k)T(L,p -l,k + 2). (4.5) 

This proposition will be proved in Section 

The recurrence (|4.5p is equivalent to the equation with which we began this paper: (|4.5p 
becomes the discrete Hirota equation (jl.ip or octahedron recurrence, when we change 
variables to n = L — p — k, i = 2p + k — L, j = p + k, and T(L,p, k) = f(n, i,j): 



f(n,i,j)f(n - 2,i,j) = f(n- l,i - 1, j)f(n -l,i + 

+ T 2 f(n-l,i,j-l)f(n-l,i,j + l). 
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FIGURE 6. A PCSSCPP of size L = 13 with a puncture of size L — 2p = 5. 
PCSSCPPs are enumerated by non-intersecting lattice paths, as indicated 
by the red lines. The paths in this picture define a natural subset of all 
PCSSCPPs, and correspond to those in Figure EJ In the East region, 
diagonal steps below the bisecting line (green online) are assigned a weight 
t 2 . All other steps are assigned weight 1. The value of k, k = 3 in this 
figure, determines the position of the bisecting green line. 

5. Fully packed loop diagrams and the Razumov-Stroganov conjecture 

In this section we will discuss another combinatorial interpretation of the solutions of 
the qKZ equation. At r = 1 the qKZ equation is equivalent to an eigenvalue equation for 
the transfer matrix of the 0(1) loop model. In this context a relation was conjectured in 
|19j (as a variant of a similar conjecture by Razumov and Stroganov [221 fl] [23]) between 
the solutions in Section [3] for L even and refined enumerations of vertically symmetric 
alternating-sign matrices (VSASMs). For L odd there is a connection to related objects, 
see below. The qKZ equation is a generalisation of the 0(1) eigenvalue equation to t ^ 1, 
as found by Pasquier [18], and Di Francesco and Zinn- Justin [8]. 

The Razumov-Stroganov conjecture can be roughly described as follows. First a simple 
bijection of VSASMs to fully packed loop diagrams is made, see e.g. [20]. To each such 
fully packed loop diagram is associated a Dyck path; this will be described below. For a 
given Dyck path a there correspond many FPL diag rains, and. their number is precisely ip a 
at r = 1. For general r there is as yet no interpretation of ijj a as a weighted enumeration 
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of FPL diagrams or VSASMs. However, as we will show in Proposition [3] at the end of 
this section, there is such an interpretation for the total r-normalisations S-(2n, n — 1) of 
the qKZ solution. We will see that S+(2n, n — 1) and S-(2n — 1, n — 1) are also related to 
ASM generating functions. 

As shown in the previous section, the normalisation S_(2n, n — 1) is also equal to a 
generating function of r 2 -weighted PCSTCPPs with a small puncture of size 2, establishing 
a connection between such plane partitions and weighted VSASMs. A similar connection 
was considered in [7j . Other equalities between generating functions of weighted cyclically 
symmetric plane partitions and generating functions related to weighted enumerations of 
symmetric ASMs were conjectured by Robbins and Kuperberg \2A\ I13| . We will prove 
some of these in Section [8l 

Finding the r-statistic on FPL diagrams that gives rise to an interpretation of ip a for 
r / 1 should provide clues for explicit bijections between (P)CSTCPPs and symmetry 
classes of ASMs, as well as for an explicit bijection between TSSCPPs and unrestricted 
ASMs. In fact, in this section we will make another small step towards such a bijection by 
elucidating the role played by the parameter p which determines the size of the puncture 
in PCSTCPPs. 

5.1. Fully packed loop diagrams. For L even, consider the set of vertices of a piece of 
square lattice of size L/2 x L. In the case of L odd, take a piece of size (L — l)/2 x (L+ 1). 
On this set of points, draw bonds such that each internal vertex has exactly two drawn 
bonds. On the boundary we impose the condition that every other ingoing bond on the 
left hand side, bottom and right hand side is drawn, starting at the top left vertex. In 
this way, the drawn bonds form closed polygons or connect outgoing bonds to each other, 
see for example Figures [7] and [8l From Kuperberg's work [13] it follows that FPL diagrams 
of even size are equinumerous with VSASMs. Those of odd size are conjecturecH to be 
equinumerous with CSTCPPs [19] . 




L/2=6 



4 5 6 7 8 9 



L=12 

Figure 7. An FPL diagram of size L = 12 with L/2 — p = 2 loop lines 
connecting loop terminals 1,2 with L,L — 1 respectively. 

We will label the outgoing bonds by successive integers, as in Figures [7] and El As each 
outgoing bond is connected to another outgoing bond, FPL diagrams can be naturally 

2 These boundary conditions and the shape of the lattice are of particular relevance in this paper, but 
more general FPL diagrams can be considered, see e.g. [3]. 
*^To our knowledge no proof of this assertion exists. 
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L -2p -1=4 




L+l=14 



(L-l)/2=6 



Figure 8. An FPL diagram of size L = 13 with [L/2\ - p = 2 loop lines 
connecting external loop terminals 2, 3 with L, L — 1 respectively. 



labeled by link patterns, or equivalently, two-row Young tableaux or Dyck paths. For 
example, the diagram in Figure [8] has link pattern ((((())(()))) which is short hand for 
saying that 1 is connected to the top of the diagram, 2 is connected to 13, 3 to 12 etc. In 
general, to each link pattern correspond many FPL diagrams. This information can also 
be coded in two-row standard Young tableaux. The entries of the first row of the Young 
tableau correspond to the positions of opening parentheses '(' in a link pattern, and the 
entries of the second row to the positions of the closing parentheses ')'. The FPL diagram 
of Figure [8] carries as a label the standard Young tableau given in Figure [9j 



1 


2 


3 


4 


5 


8 


9 


6 


7 


10 


11 


12 


13 





Figure 9. Standard Young tableau corresponding to the FPL diagram in Figure El 

Yet another way of coding the same information is by using Dyck paths. This will be 
useful when making a connection to the results of Section [3l Each entry in the first row of 
the standard Young tableau represents an up step, while those in the second row represent 
down steps. The Dyck path corresponding to Figure [9] is given in Figure [TUl 




Figure 10. Dyck path corresponding to the FPL diagram in Figure [8] and 
the standard Young tableau in Figure [8j 



5.2. Subsets of FPL diagrams and punctured plane partitions. As discussed above, 
FPL diagrams of even size are equinumerous with VSASMs, and it follows from the results 
of Section HI that they are also equinumerous to PCSTCPPs with a small puncture of size 
2. Here we shall formulate a refined correspondence between subsets of FPL diagrams and 
PCSTCPPs. 
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Recall the subset £>l,p of Dyck paths whose local minima lie on or above height p where 
p = \_{L — l)/2j — p, see Sectional Each path in T>L tP is a label for FPL diagrams whose 
loop terminals 1, . . .p (in the case of L even) are connected to terminals L — p + 1, . . . , L 
respectively, see Figured In the case of L odd, the Dyck paths in T>l >p label FPL diagrams 
with loop connections between terminals 2, . . . ,p + 1 and L — p+l,...,L respectively, see 
Figure M 

Definition 3. A p-restricted FPL diagram is an FPL diagram whose corresponding 
Dyck path belongs to T>L tP . 

In other words, in a p-restricted FPL diagram of even size, the first L/2—p loop terminals 
are connected to the last L/2 — p loop terminals. This classification of FPL diagrams 
allows us to formulate a refined correspondence between FPL diagrams and symmetric 
plane partitions: 

Conjecture 2. The total number of p-restricted FPL diagrams of size L = 2n is equal to 
the total number S±(2n,p)\ T=1 of PCSTCPPs having sides of length 2{L — p) and 2(p + l) 
with a triangular puncture of size 2{L — 2p — 1), see Figure^ 

As a side remark we note that S±(L,p)\ T=1 at r = 1 factorises (this was proved in [T5j [3] 
in the context of PCSTCPPs) and takes the form 

c (T M _ 9(P+ 1)(L- P ) TT r(L - j + l)r((2L + 2j + 3)/6)r((L - 2j + 3)/3) 

*±^,p)\ T =i-^ JUL r(L-2j + i)r(j + i/2)r((2L-j + 3)/3) ' [ ' 

This expression is equivalent to earlier conjectured expressions for the total number of p- 
restricted FPL diagrams in the context of the 0(1) loop model and the Razumov-Stroganov 
conjecture pTT], [2Tj . 

One may ask if a generalisation of Conjecture [2] holds for t / 1. As we have seen 
in Section 14.11 the parameter r 2 is a natural weight for (punctured) symmetric plane 
partitions. One can therefore ask whether it also describes a natural statistic on ASMs or 
FPL diagrams. In the case of p = \_(L — l)/2j (p = 0) this question can be answered in 
the affirmative. It was already observed by Robbins and Kuperberg \2i\ [T3] that some r 2 - 
generating functions for ordinary CSTCPPs, i.e. without puncture, are the same as certain 
polynomials arising in r 2 -enumerations of symmetric ASMs where each —1 is assigned a 
weight t 2 . In the case of FPL diagrams this amounts to giving a weight r to every two 
consecutive vertical or horizontal steps (such consecutive steps correspond to either a +1 
or —1 in ASM language). Our expressions for the partial sums in Table [2] can be directly 
compared with Kuperberg's Table 4. Using the notations of [13], see also Section [8l we 
collect this observation with two others in the following intriguing proposition: 

Proposition 3. 

r(2n,n-l,2) =S-(2n,n-l) = A Y (2n + 1; r 2 ), 
T(2n,n- 1,1) = S + (2n,n-l) = A^ p (4n + 2;r 2 ), 
T(2n- l,n- 1,1) = r n ~ 1 5_(2n - 1, n - 1) = i$j(4n; r 2 ). 
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This proposition will be proved in Section El 

It is an open problem to find a generalisation of Conjecture [2] to arbitrary weight r, 
or a generalisation of Proposition [3] to arbitrary p, i.e. to find a correspondence between 
r 2 -weighted p-restricted FPL diagrams and r 2 -enumerations of PCSTCPPs. 

6. Proof of the determinant formula for partial sums 

In this section we prove Proposition Q] using the formalism developed in [91 [32] and 
which consists of writing integral formulae for solutions of the qKZ equation. For the sake 
of simplicity, we shall work out separately the two possible parities of the size L. 

6.1. Even size. Assume L = 2n. The following set of integrals was introduced in [9]: 

1pa u ...,a n = [1+Xi- Xj][l - Xi ~ Xj] I ■ • • I JJ 

m<n [yt - y m ][i + vi- y m ][ye + Vm] ilK^^J 1 

~t~ yi ~l~ 2/mJ , , 

= = 1 y — — = (,0-1) 

nLi n i= i [ye + Xi] UiU [ye - x i] Ui= ae +i t 1 + yt - x i] 

where the contour integrals surround the poles at Xi — 1. Here a±, . . . ,a n form a non- 
decreasing sequence of integers between 1 and L — 1. 

The relation to the solution of the qKZ system (|2.3p is as follows: up to normalization 
by a symmetric factor of the parameters x^, the Vai ... o„ are linear combinations of the 
components of the solution: 

i> ai ,...,a n =^2C ait ... >ania ij) a (6.2) 

a 

with coefficients C aij ... |0n;a that are described explicitly in appendix A of [32], and which 
we shall define here by recurrence. First, C§.§ = 1. Next, for a pair (a\, . . . , a n ; a) of 
length L, consider any local maximum i of the path a, and the new path a' obtained by 
removing the two steps before and after i. Call k the number of I such that ai = i. If 
k = 0, C a , .,a n -,a — 0. If k > 0, consider the new sequence a'-y, . . . ,a! n _i obtained from 
ai,...,a n by removing one "i" and replacing each remaining ae with: itself if < i; 
i — 1 if an = i; at — 2 if ag > i. Then C olr .. ian;a = [fc]C a / ;a / (and this definition is 

independent of the choice of local maximum). 

Since we are interested in the values of the ip a at xi = 0, let us set X{ = in (|6.1h and 
perform the change of variables ui = [1 + y n+ i-i\/[y n+ i-i\. At this stage it is convenient 
to reindex the integers as bp = L — a n+1 _£, I = 1, . . . ,n, and to define ^f, li ... i b„ := ip ai ,...,a n , 
so that 



^ax,...,a n = VV...,fc„ =*■■■*] 



, 27riu^ f 



K£<m<n 



X 



(1 - l^lt„ 

(6.3) 



(u m - u^)(l + rn m + u e u rn )(T + u t + u. 

l<£<m<n 
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Here the normalization is chosen in such a way that ipi,2,...,n = ?Pn = r n(n_1)/2 (note that 
in this case the integrals are trivial and can be performed by simply setting ui = in the 
numerator of the integrand). 

We now consider specific 'tjjb 1 ,...,b n which will reproduce our sums S±(L,p). Fix a non- 
negative integer p, and let p = n — 1 — p. Consider sequences (&i, . . . , b n ) of the form 

U i<t<p+i 

1 ~ \2£ -p- 1 - tt-y-x p+2<£<n 
where ei,...,e p €{0,l}. We have the following 
Lemma 1. 

$bi,...,b n = '01,...,p+l,p+3-ei,...,2n-p-l-e p = ^ Ipa 

V^,Qi_p_2«-l<ai-p-2< iff Q=l 

Proof. We shall proceed by induction. Fix a sequence of integers be as in the lemma, that 
is in terms of the mirror-symmetric sequence a,£, 



2£ + p - 1 + e p+1 _ e \<£<p 
i+n-l p+l<£<n 



Note that an > p + 1. Let a be a Dyck path. Consider a local maximum i of a. One can 
always assume i < n. There are two cases: 

1. i < p. In this case, a ^ T^L,p- We find immediately that there are zero ai = i, so 
that the coefficient is zero. 

2. p<i<n<L — p. Call £ = [(i — p + l)/2j . There are four cases depending 
on the parity of i and the value of e p +i_^. If e p+ i_^ = and i = 2£ + p ^ ag, there 
are no ag equal to i so that the coefficient is zero. Since i is a local maximum, we 
have a£_p„ 2 (p+i-f)-i=i-i < a L-p-2(p+i-£)=i satisfying the inequality in the summation 
of lemma 1 (despite e p +i-e = 0). Similarly, if e p+ i-e = 1 and i = 2£ + p— 1 ^ a£, there are 
no a£ equal to % so that the coefficient is zero, and a L _p_ 2 ( P +i-e)-i=i > aL-p-2(p+i~e)=i+i 
violating the inequality in the summation of lemma 1 (despite e p+ i_^ = 1). In the other 
two cases, we have = i and the equivalence in the summation of lemma 1 is valid. We 
can then apply the definition by recurrence of the coefficient C aii ... jan;Q ; the new sequence 
a'y,... ,a' n _ 1 is exactly the same type as a\, . . . ,a n , that is defined by the same ei with 
e p+ i_^ skipped. On the other hand it is clear that the other conditions on a in the 
summation of lemma 1 are equivalent to the conditions on a' (a with the local maximum 
removed) with the new sequence a'i, . . . , a' n _i. One then uses the induction hypothesis to 
conclude. □ 

Example 1. Consider the two sequences with L = 6, p = 1 (p = 1). We find 

$l,2A = i> , 2 =2r 2 (l+r 2 ) 

^1,2,3 =tp /\ , = T 3 
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In size L = 8, with p = 2 (p = 1), 

= 4, 



6r 3 + 21r 5 + 18r 7 + 5r 9 



5r 4 + 7r b + 3r s 




3r 5 + 3r 7 



XTie edges in red are those whose labels (counted from right to left) appear in the sequence 
of integers iff the edge is a down step ( in fact the first and last p + 1 edges are fixed by the 
fact that the paths are in T>l iP ). 

Note that taken together, the various sequences for a given L and p reproduce the full 
set of paths of T>l >p - Furthermore, by direct computation using formula (13. ip (grouping 
together pairs ctL-p-2i-i and ai-p-2i in the sum and using_aL-p-2j-aL-p-2i-i = 2ej-l), 
it is easy to check that all ip a that contribute to a given ■0i,...,p+i,p+3-ei,...,2n-p-i-e p have 
the same integer c a:P = Y2i=i e i- 

We thus define 

S(L,p\t)= ^ ^= ie ^l,..,p+l,p+3-e 1 ,..„2n-p-l- ep (6.4) 
ei,...,£ p e{0,l} 

and claim that S±(L,p) = 5(L,p|r ±:L ). 

Using (|6.3p . we now obtain the following integral representation for S(L,p\t): 



S(L,p\t) = f---f(j[ 



n 



due(l + tue) 



U-p-l 



(1 - u t u m ) 



\<Km<n 



(6.5) 



X (u m - U e )(l + TU m + U £ U m )(T + U£ + U m ) 

l<£<m<n 

The first integrals over u±, . . . , ttp+i can be performed successively by simply setting the 
corresponding variables to zero. The result, after shifting the indices of the variables, is: 

' du e (l +tu e )(l + Tu e )P +1 {T + u e )P +1 
=i 



s(L, P \t)= /•••/n 



2m uf 



(1 - u t u r 

l<l<m<p 



X (u m - U £ )(l+TU m + U e U m )(T + Ui+U m ) 

l<£<m<p 

Next we use the following lemma: 



(6.6) 
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Lemma 2. If AS designates antisymmetrisation: AS(f(ui, . . . , u p )) = ^o-es ( — l) CT /( n o-(i)> ■ 
and (■••)<() means keeping only non-positive powers of a Laurent polynomial in the vari- 
ables up, then the following equality holds: 



(1 - upu m ) AS Yl u e 2£+1 U + u * u ™ + TU r 



^l<l<m<p 



l<i<m<p 



" 



Asjn^v+^Y- 1 ) =n«r 1 n (^-^xt+v'+o (6.7) 

/ f=l !<^<m<p 



^=1 



This is equivalent to formula (4.5) of [32J. It is a slightly stronger version of the propo- 
sition of [30], and can be proved along the same lines. We use it to symmetrize the 
integrand: 



S(L,p\t) 



1 



yr du e {l + tup)(l + rupf +1 (T + upf +1 

(u m - ui)(t + ue + u m )(u^ - uj 1 )^ + UJ 1 + u~}) (6i 

!<£<m<p 



Noting that T[i<p <m < p (u m ~ u t) (t + up + u m ) is just the Vandermonde determinant of the 
u/(t + up) and similarly for the other factors, we can finally pull the determinants out of 
the integral, resulting in: 



S(L,p\t) = det 

l<e,m<p 



du 
2iriu 



(1 + tu)U l - m+P {T + U) £+ P(T + U" 1 ) 



-1 \m+p 



(6.9) 



For general t, by using the binomial formula we can evaluate this to be 



S(L,p\t)= det 

l<£,ni<p 



p+2m+2€-2r-l ( l + P\ (J m +P 



r — i ) V V 2m — r 



m + p 
2m — r — 1 



where we recall that p = L/2 — 1 — p. 

At t = r this expression simplifies slightly: 



S+(L,p) = det 

l<t,rn<p 



V- T p+2m+2£-2r f£+p\fm+p+l\ 

^ \r - i) V 2m - r J 



as well as at t = r : 



S-(L.p) = det 

l<£,m<p 



.p+2m+2£-2r 



r — ^ / V 2?ra — r 



(6.10) 



(6.11) 



(6.12) 



The summation over r is such that only a finite number of terms is non-zero; in practice 
a possible range is < r < 2p. One can check that formulae (|6.11l6.12p match the 
expressions given in proposition Q] for L even. 
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6.2. Odd size. Assume L = 2n + 1. The reasoning being exactly identical to the case L 
even, we only provide the key formulae. The starting point is formally the same integral 
formula as previously: 



\\ [I + Xi - Xj][l - Xi - Xj] j> ■■■ j> \\ 



log q dy e 
q — q" 1 7ri 



i<i<j<L " " e=i 

Ul<t<m<nlyi ~ + Vt- Vra][yt + Vm] X[l<l< m <r\ l + Vl + V™} 



(6.13) 



but with an odd number of parameters Xi, so that it produces a slightly different expression 
when the Xj £1X6 set to zero: 



' l Pai,...,an — ^6i,...,6„ 



n 



dn^ 



(1 - u^u m )(l + TU m + u^« m ) 



KKro<ti 



X 



(«m - + + Urn) 



l<£<m<n 

Define, for p a non-negative integer, p = n — p and 



(6.14) 



S f (-E'.P|*)= J] * E ^ iei, 01,...,p,p+2-e lv ..,2n-p- ep 
ei,...,e p e{0,l} 



(6.15) 



so that S±(L,p) = S(L,p\t ±1 ). We obtain the following integral expression for S(L,p\t): 



dui(l + tue) 



2iriu 



2l-v 



(1 - u e u ri 

l<e<m<n 



X JJ(1 + TU£ + uf) (U m ~ Ui)(l + TU m + U e U m ){T + U£ + U Tl 

1=1 l<£<m<n 

We integrate over the first p variables and reindex the remaining ones: 



(6.16) 



s(L, P \t)= /•••/n 



du e (l+tUl){l+Tuif(T + Uif 



2ni uf 



(1 - u £ u r 

l<e<m<p 



+ TU £ + u|) (ti m - U t ){\ + TU rn + U e U m )(T + U t + U Tl 

t=\ l<e<m<p 



(6.17) 
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We use lemma 2 and pull the determinants out of the integral as before: 



S(L,p\t) 



1 



n 



du t (l + tu e )(l + ruif{r + utf{\ + rue + uj) 
27ri v?, 



det 

l<£,m<p 

det 

l<l,m<p 



(u m - U£)(t + u £ + u m ){u rr } - u t 1 )(t + u e 1 + u n }) 

l<£<m<p 

du 



2iriu 
du 

27T1U 



(1 + tu)(l + TUg + 



l\m+p— 1 



(i + 



-m+p / 



(t + ^+p-^t + u- 1 ) 



hm+p 



(6.18) 



In the last line we wrote (1 + ru + u 2 ) /u = (t + u" 1 ) + u and noted that the second term 
reproduces the column m — 1 of the matrix and thus can be subtracted without changing 
the determinant. 
One finally obtains 



S(L,p\t) = det 

i<e,m< P 



T P+2m+2£-2r-3 



l+p-1 

r-l 



m + p 
2m — r 



+ t 



m + p 
2m — r — 1 



where we recall that p = (L — l)/2 — p. 
At t = r this simplifies to 



S + (L,p) = det 

l<£,7Tl<p 



E 



p+2m+2f-2r-3 



£ + l\/m + p+ l 



r-4 



2m — r 



whereas at t = r 



-l. 



S-{L,p) = det 

l<^,m<p 



E 



,p+2m+2^-2r-l (£ + P\ f TTl + p 



(6.19) 



(6.20) 



(6.21) 



K r — £ J \2m — r , 
thus reproducing the expressions of proposition Q] for L odd. 

7. Proof of the bilinear recurrence relations for T(L,p,k) 

In this section we prove Proposition [2j To simplify the formulae of this section we 
introduce k' = L — 2p — k, so that definition (|4.ip reads 



T(L,p,k) 



det T lm 

l<£,m<p 



T em (k,k')= 



t + k-l\fm + k' 

2m — i — r 



'2r 



(7.1) 



Here the limits of summation in r are automatically fixed by the conditions ( n ) = 
Vr < and Vr > n. We begin with a derivation of yet another determinant formula for 
T(L,p,k). 



Lemma 3. 

where 



T(L,p,k) = det Utm, (7.2) 
l<z,m<p+± 

Ui,i(p)= (-l)'- 1 ^ 1 -*), Ue, m+1 (k,k')=T em (k,k'-l). (7.3) 
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Proof. To check the identity det U = det T we shall perform linear transformations of the 
matrix U not affecting its determinant. First, we combine adjacent rows of U with the aim 
to set to zero all components of the first column, except the last element Up+n = (— 

,m+l • 

Then, by decomposing the determinant of the resulting matrix along the first column we 
find 

det Uf m = det Vi m . 

l<£,m<p+l l<£,m<p 

The lemma now follows by noticing that the rows of T are linear combinations of those 
of V, 



l + k-l\fm + k'-l\ ( £ + k \fm + k'-l 



^ = E L./J r + 



, r 

2m — i — r J V T I \2m — I — r J V r — 1 



2r 



T 2r 



~ ^\\2m- £-r)\ r ) + \2m - £ - r - l) \ r-1 

= T tm {k, k') + T 2 T^ m ^i(k, k'), 

where we have used Pascal's rule (?) = ( n ~ 1 ) + (^J) to go from the second to the third 
line. Hence we find 

m— 1 

^ (— r 2 ) J Vg, m -j = Tim, 

from which we conclude that det V = det T. □ 



Our derivation of formula (j4.5j) is based on the use of the following particular example 
of a Pliicker relation for determinants (for the general case see [27] ). Consider a pair of 
n X n matrices A^ m and B^ m . Denote by Ag the £-th row of the matrix A and introduce 
notation 

detA= \A\, A = [Ax,..., An], 
In this notation the Pliicker relation reads 

n 

\A\ \B\ = l^i) ■ ■ ■ j Ai-ii Bj\ x ... , Bj-i, A n , Bj + i, . . . ,B n \, (7.4) 
i=i 

where the sum is taken over permutations of the last row of A with each row of B. 

We now take n = p + 1 and, recalling the definition of the matrix U in (jT.3|) . substitute 
for t4 and -B the following matrices: 

A(k, k') = U(k, k') = [Ui, ... , U p+1 ] , (7.5) 
B(k,k') = [S u U 2 ,...,U p ,S p+1 ], (7.6) 

where 8i m is Kronecker's delta. 

By Lemma [3] we have \A\ = T{L,p,k). To calculate the determinant of B we expand 
along its top and bottom rows and then notice that the first column of the resulting 
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(p — 1) X (p — 1) matrix is also of <5-type: 

^, 2 (A:,fcO = T,, 1 (A ; ,^-l) = ^Q+^^(^r 2 '- = ( 5,, 2 for i > 2. 

So we calculate 

|.B| = det Ue +2 .m+2(k,k') = det T^ +2 , m +i 1) 

l<«,m<p-2 ' l<£,m<p-2 

det 2i m (fc + 2, fc') = T(L - 2,p - 2, fc + 2). 

l<£,m<p-2 

Now let us consider the right hand side of the relation (17. 4h . Here only permutations 
with the top and the bottom rows of the matrix B give non-vanishing contributions. 
Permuting the last row in A and the first row in B gives the following factors: 

\U u ...U p ,Si\ = (-If det U e>m+1 (k,k') 

l<£,m<p 

= (-l) p det T em (k,k' -I) = (-l)PT(L-l,p,k), 

l<£,m<p 

\U p +i, U 2 , ■ ■ ■ U p , 5 p+ i\ = (— l) p ~ 1 |{7 2 , • • • > £^>j t^p+i| 

= (-If det U im {k + 2,k' + 1) = (-l) p T(L- l,p- 2, k + 2), 

l<£,rn<p—l 

where in the last calculation when passing to the second line one i) expands the deter- 
minant along the first column noticing that E/f+i 2 = and ii) redefines matrix indices 
using the identity Ui + 2,m+i(k, k') = Ui m (k +2, k' + 1) V m > 1. 

Permuting the last row in A and the last row in B gives the following factors: 

\U 1 ,...,U p ,6 p+ i\ = \U U ...,U P \ = t 2 det U tm {k,k') = r 2 T(L-2,p-l,k). 

l<£,m<p 

The factor r 2 is extracted from the p-dependent first column of the matrix U : Ue : i(p+1) = 
r 2 U^{p); 

\5i,U 2 ,- ■ -U p+ i\ = det Ui + i, m+ i(k,k') 

l<£,m<p 

= det U l+2 , m+2 {k,k') = det T lm {k + 2,k') = T{L,p - l,k + 2) . 

l<£,m<p-l ' l<£,m<p-l 

Here when passing to the second line of the calculation we first expand det C/^ + i im+ i along 
the first column C^+i,i = and then redefines indices of the matrix T: I) +2im+1 (/c, k! — 
l)=T tm {k + 2,k'). 

Thus, the Pliicker relation (|7.4j) for the matrices A and B defined in (17. 5p and (|7.6|) 
produces the equality (14. 5p . 



8. Proof of t -enumeration of ASMs 

In this section, we prove Proposition [3l We start from the determinant formulae of [13] 
for the enumeration of various symmetry classes of Alternating Sign Matrices, and reduce 
them to our own determinant formulae for S±. In what follows, we keep Kuperberg's 
notations (even though they are non-standard), to ease the comparison of formulae with 
|13j . In each case, one starts with configurations of the six- vertex model, which for certain 
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particular boundary conditions are identified with Alternating Sign Matrices in various 
symmetry classes. There are three distinct Boltzmann weights for the six-vertex model, 
taking into account Z2 symmetry, and they are parametrised as w a = ax~ l y — a~ 1 xy~ 1 , 
uib = axy~ x — a~ 1 x~ 1 y, w c = a 2 — a~ 2 , x and y being row/column spectral parameter and 
a a global parameter. In the ASM language, a weight w c is assigned to a ±1, and weights 
w\, and w & are assigned to zeroes. In the end we must take the homogeneous limit where 
w & = u>b and w c /w & = r = —q — q : this ensures that adding a —1 to an ASM, that is 
two extra vertices of type c (adding a —1 also increases the number of +1 by 1), produces 
a weight r 2 . This is achieved by setting all spectral parameters to 1 and a = —q. Similar 
parameters, called b and c, which are related to boundary weights will be used below. 

To be self-contained, we will give the matrices of Kuperberg relevant to this paper. It 
is useful to first define the functions o and a by 

o~(x) = x — x~ , a(x) = a (ax) a (a /x). 

The relevant matrices then are 

Mu(n; x, y\j = _,},__ , - ,} T , (8.1) 



Muu(n;x,y)ij 



a(xi/yj) a(xiyj) 
a(b/yj)a(cxi) a(b/yj)cr(c/xi) a(byj)a(cxi) 



a(axi/yj) a(a/xiyj) a(axiyj) 

+ crjby^o-jc/xi) 
a(ayj/xi) 

M+ T (n; x, y)ij = 1 + 1 ■ (8.3) 
a{ayj/xi) a{axi/yj) 

Since the equalities of Proposition [3] are known to be true at r = 1 (they are conse- 
quences of the various relations between enumerations of ASMs found in |13| . as well the 
relations between ASMs and PPs discussed in [9]), and since both sides are easily checked 
to be polynomials in r of the same degree, we can safely drop various trivial factors in the 
calculation, keeping only the determinant itself as well as factors that become singular in 
the homogeneous limit. 



8.1. First formula. In |13J . it is explained how VSASMs are a special case of UASMs; 
more precisely, the partition function of VSASMs can be obtained from the more general 
one of UASMs by tuning a certain boundary parameter. It is also noted there that the 
boundary parameter only enters the formula for the partition function in prefactors, and 
not in the determinant itself. In particular the enumeration of VSASMs and UASMs are 
essentially the same. We shall therefore write directly the partition function for UASMs 
of size 2n without all the regular prefactors; using (|8.ip it takes the form: 

-ZuASM OC 



A(x 2 )A(y 2 )A*(x 2 )A*(y 2 ) 

1 



det 



i<M<« I (a 2 xf - y 2 )(a 2 y 2 - xf)(a 2 - x 2 y 2 )(l - a 2 x 2 y 
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where A stands for the Vandermonde determinant, e.g. A(x 2 ) = r[i<j( x j ~~ x a)'i an d 
A*(x 2 ) = n-i<j (1 ~~ x i x J)- The X{, t/j are spectral parameters which will eventually be set 
to one. 

We use the following integral representation: 
1 



(a 2 x 2 — y 2 )(a 2 y 2 — x 2 )(a 2 — x 2 y 2 ){\ — a 2 x 2 y 2 ) 

a 3 (a 2 - l) 2 f du ^(n + a + a" 1 ) 



-, (8.4) 



(a 2 + 1) x 2 y 4 (x 2 — a 2 )(a 2 x 2 — 1) J 2it\ fi(u, x)fi(u, 1/ x) fi(u, a/ y) n(u, ay) 
where 

x) = a(a + u) — (1 + au)x 2 . (8-5) 
The contour of integration in (|8.4p surrounds the x-dependent poles but not the y- 



dependent ones. This identity can be checked directly by residues. The various prefactors, 
as well as the integral sign, can be pulled out of the determinant and we thus find 

. dui(ui + a + a" 1 )^ / du n (u n + a + a -1 )u n , 

Ajasm oc 



A(x 2 )A(y 2 )A*(x 2 )A*(y 2 ) J 2vri J 2vri 



det | — A — — - ) det 



^fi(ui,Xj)n(ui,l/xj) J \n(ui,a/yj)n(ui,ayj) 
In order to compute these determinants, we perform the following change of variables: 
x _ (1 ~ x 2 ) 2 y _ (1-y 2 ) 2 



(1 - a 2 x 2 )(l - a~ 2 x 2 ) (1 - a 2 y 2 )(l - a~ 2 y 2 ) 

and use the factorizations 

fi(u, x)/j,(u, 1/x) = (a 2 — x 2 )(a 2 — x _2 )(l — Xu(a + a" 1 + u)), 
fi(u, ay)n(u, a/y) = a 2 u 2 (a 2 — y 2 )(a 2 — y~ 2 )(l — Yu~ 1 {a + o _1 + u~ 1 )). 

Note that a + a" 1 = r. Again one can get rid of the trivial factors and obtain 

1 fduiiux+r) fdu n (u n + T) 

^UASM oc 



A(x 2 )A(y 2 )A*(x 2 )A*(y 2 ) / 27riiti J 2m u n 



det I r— — — ) det 



1 - Ui {r + m)Xj J V 1 - u- % (t + u^ 1 )Yj J ' 

where the contours of integration surround the Xj-dependent poles. The determinants are 
now of Cauchy type and can be evaluated exactly: 

A(X)A(Y) /duiOi + r) Jdu n (u n + T 

^UASM OC 



A(x 2 )A(y 2 )A*(x 2 )A*(y 2 ) J 27riui J 2vriu 

A(«(l + ™))A(«- 1 (1 + ru- 1 )) J] „,_,.,^ N ,| 



1 (i - ui(r + - «r + «r )*;•) 
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The Vandermonde determinants outside the integral cancel each other, leaving only a 
regular part, due to 

(l-a 2 ) 2 (x?-xf)(l-x 2 x 2 ) 

Ai - A; 



J a 2 (l - a 2 x 2 )(l - a" 2 x 2 )(l - a 2 x 2 )(l - a" 2 x 2 ) ' 

At this stage one can take the homogeneous limit, that is set Xj = y% = 1, or X$ = Yi = 0, 
which results in the simple expression 

9 . f dui(r + Ui) f du n (T + u n ) . , , . . . . i . , ,. 

A v (2n + l;r 2 ) oc i — 1J - ■ ■ ■ <b — ^ —A(u(t + u))A(u (t + u )), 

J 2mm J 2niu n 

where the integrals are performed around zero. We can once again exchange determinant 
and integral sign and find 

A v (2n + l;r 2 ) oc det / -^ u ^ m {r + u) 1+1 {t + U" 1 )™ 

0<e,m<n-l J 2lT\U 

The first column of this matrix being (1,0,...), we can restrict the range of the indices to 

A v (2n + l;r 2 ) oc det / -^ u ^ m { T + n) m (r + u" 1 )™ 
l<£,m<n-l J 2lT\U 

which is identical to (|6.9I) with t = r" 1 , p = n — 1, p = 0. 

8.2. Second formula. The reasoning is exactly the same; only the matrix elements of 

(2) 

the determinant are slightly modified. In the case of ^vhpasM' we mus t use t ne matrix 
Mtjtj defined in (18. 2D with parameters b = 1/a, c = a. In this case we use the following 
identity to represent the matrix elements in factorised form: 



a 2 y 2 (l + x 4 ) + ((a 4 + a 2 + l) (1 + y 4 ) - (a 2 + l) 2 y 2 ) x 



„2 

„2 



xy(a 2 x 2 — y 2 )(a 2 y 2 — x 2 )(a 2 — x 2 y 2 )(l — a 2 x 2 y 2 ) 

= (l-a 2 ) 



2 . 2 l u(l + ti(a + a x )) 



27ri x)(ti, 1/x)/j,(u, a/y)fj,(u, ay) 

(2) 

We thus find the following expression for the "partition function" -ZyHPASM (which is 
really a ratio of the partition function of UUASMs by the partition function of UASMs) 

(2) A(X)A{Y) r dujXl + rui) f du n (l + ru n ) 

VHPASM a A(x 2 )A(y 2 )A*(x 2 )A*(y 2 ) / 27rim "'/ 2vrin n 

1 



i (1 - Ui (r + u i )X i )(l - nr!(r + u^Yj) ' 



A(«(l + ™))A(n~ 1 (l + tu' 1 )) \\ - 
and in the homogeneous limit, 

4 2 HPASM(2n + 1; r 2 ) oc f -f ^^A(u(r + „))A(™^(r + O) 

oc det /^^- m (r + ^(r + «- 1 r +1 . 

0<£,m<n-l 7 27T1 

Once again one can remove the first line and column, and we recover (|6.9p with t = r, 
p = n — 1, p = 0. 
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8.3. Third formula. The generating function A^(n; t 2 ) is denned as one of the factors 
of A^(2n;r 2 , 1), the other being n;r 2 ,l,l). To compute this generating function 
we need Z^ 2 \2n; (x, x _1 ), (y, y -1 )) which is defined in terms of the matrix Mjtp, see 
^3|) . The matrix Mg T (2n; (x, x _1 ), (y, y -1 )) commutes with 



P 



/„ 
In 



and thus can be brought to block-diagonal form with blocks 

M+^ ± (n; x, y)ij = M+ T (ra; x, y)ij ± ^h T (™; iT 1 )^- 

These matrices are readily identified with the matrix Mtjtj with b = c = i and b = c = 1 
respectively: 

M^f(n; x, y)ij = 1 M UT j(n; x, (6 = c = i), 

cr(ix)<7(iy) 

M^ + (n; x, y)ij = - , . Muu(n.; x, y)y (6 = c = 1). 
The further reasoning is again analogous to that in Section 18.11 and in order to compute 

(2) -» 

the homogeneous limit of det M^(n; if, y) we now use the identity 
(a 2 + l) 2 x 2 y 2 - a 2 (x 2 + y 2 )(l + x 2 y 2 ) 2 



xy(a 2 x 2 — y 2 )(a 2 y 2 — x 2 )(a 2 — x 2 y 2 )(l — a 2 x 2 y 2 ) 



a 



22 , du 



(l-a^r* -, T7 r-— r. (8.6) 

Z7ri x)(ii, l/x)fi{u,a/y)fj,[u, ay) 



Thus we find 



- (2 ) ApQAflQ /■ dm /■ du 

Ztttt oc 



uu A(x 2 )A(y 2 )A*(x 2 )A*(y 2 ) J 27riui J 2mu n 

A(«(l + ™))A(n" 1 (l + rn" 1 )) J] * - V , - hvV 

and 

ig) (4n; r 2 )<xf^--..f ^A(u(r + u))A(^(t + n" 1 )) 
oc det /— /-"(r + nWr + n-'r. 

0<£,m<n-l y 27Ti 

Removing the first line and column, we recover (|6.18p with t = r _1 , p = n — 1, p = 0. 
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